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This paper is concerned with statistical properties of a gas of gp-bosons 
without interaction. Some thermodynamical functions for such a system in 
D dimensions are derived. Bose-Einstein condensation is discussed in terms 
of the parameters q and p. Finally, the second-order correlation function of 
a gas of photons is calculated. 



1. In the last twenty years, different types of deformed bosons have 
been introduced [1-8]. There use for obtaining realisations of quantum 
algebras is now widely widespread (e.g., see [4-6,8]). We shall deal here 
with gp-bosons [6-8] associated to the annihilation (a), creation (a + ) and 
number operator (N) characterized through their action 

a \ n ) = \l[[ n ]] qP \ n ~ !>> a+ \ n ) = \J[[ n + 1]W l n + !>> N\n)=n\n) (1) 

on the ordinary Fock space T := {|n} : n G N}. (The operators a, a + and 
TV span a gp-deformed Weyl-Heisenberg algebra.) The gp-deformed number 
in (1) are defined by 

[WW := V" < 2 > 

— P 

The particular cases p = 1 and p = g _1 correspond to g-deformations 
mainly used in the mathematical and physical literature, respectively. In 
the general case, Hermitean conjugation requirements demand that (q,p) G 
R x R, or (q,p) G C x C with p = q, or (g,p) G S 1 x 5 1 with p = q = q~ x 
[9], 

Numerous works have been devoted to the study of the statistical ther- 
modynamical properties of a gas of deformed bosons (mainly g-bosons) 
without interaction [10-21]. (When writing this Letter, we became aware 
of a paper [21] dealing with gp-bosons.) In particular, the phenomenon of 
Bose-Einstein condensation of deformed bosons has been studied by several 
authors [10,13,14,16,20,21]. Furthermore, the deformation of the correlation 
function of order two for a gas of photons has been recently investigated by 
Man'ko et al. [15]. Most of the preceding works concern q- and gp-bosons in 
D = 2 or 3 dimensions and are restricted to q G R [10-20] and (q,p) G RxR 
[21]. 

We derive in this Letter the (thermodynamical) properties of a gas of 
gp-bosons in D dimensions for parameters q and p varying in five domains. 



The Bose-Einstein condensation of qp-bosons is discussed in a new approach. 
Finally, we obtain the correlation function for a gas of photons described 
by gp-bosons. 

2. We start with the Hamiltonian (in a second-quantized form) 

H:=Y,H k , H k := (E k - i2)N k , (3) 
k 

which describes a gas of gp-bosons without interaction. In eq. (3), \i is the 
chemical potential while E k and N k are the kinetic energy of a gp-boson 
and the number operator for the gp-bosons, in the k mode, respectively. 
The gp-analogue of the Bose factor for the k mode is then 

(/ fe W:=|tr(e-^a+a fc ), (4) 

where Z := tr (e _/3H ) is the partition function and (3 = (/csT) -1 the recip- 
rocal temperature. As a trivial result, we have 

Z = UT^ r,:=j3(E k -fj). (5) 
k 1 e 

The partition function Z is thus independent of the deformation parameters 
q and p. Then, the Bose factor (fk) qP reads 

(fk) gP = (1 - e-v)tr(e-^[[N k }} qp ), (6a) 

which can be shown to converge in each of the following five cases : (i) 
(q,p) G R+xR+ withO < q < e _/3 ^ and <p < e~^; (ii) (q,p) G R + xR+ 
with p = q~~ l and e^ M < q < e~^ M ; (hi) (q,p) G R + x R + with p = 1 and 
< q < e _/3 ^; (iv) (q,p) G C x C with p = q and < \q\ < e _/3M ; (v) 
(q,p) G S 1 x S 1 with p = q = q~ x . 



At this stage, it is important to justify the choices made for eqs. (3) 
and (4). The gp-bosons are not new particles. They are still bosons. Thus, 
(3) and (4) correspond to the basic formulas for the statistical mechanics 
of bosons. Here, we do not deform the basic formulas. The trace and the 
exponential in eq. (4) are nondeformed functions. We only use (through 
the formalism of second quantization) creation and annihilation operators 
that satisfy #p-deformed commutation relations. The form chosen for 
in eq. (3) is nothing but the usual one where the number operator takes its 
eigenvalues in N. The deformation cannot enter the theory at this level by 
replacing Nk by a^a,k in (3) since we would obtain an "occupation number" 
a t a k = [[Nk]]qp that would not have its eigenvalues in N. On the contrary, 
we must keep a\ak rather than Nk in (4). Indeed, should we use Nk in 
place of a^ak in (4), we would get a nondeformed value of the Bose factor. 
As a conclusion, eqs. (3) and (4) constitute a reasonable starting point for 
a (minimal) deformation of the Bose distribution. 

Following ref. [9], we can use the parametrisation 

q = e ^cosT e +^sinr^ ^ _ Q <p cos T Q -<p sin r ^ 

for (q,p) G R x R and the parametrisation 

q = cos T e +i ^ sin T , p = e^ cos T e~ i ' p sin T (8) 

for (q,p) E C x C or (q,p) G S 1 x 5 1 , where ^ER* and -tt/2 < r < tt/2, 
so that the limiting situation q = p — 1 corresponds to cp — > 0. 
The obtained expression for (fk)qp is 

(t \ 9-11 p - 1 1 

(fk) qP = + (9a 

q — pe 7 ! — q p — q e^ — p 

or alternatively 

e v - 1 

(/fe)w = ( e v - q)(ei - p) 1 (10) 



which clearly exhibits the q <-> p symmetry. In the (five) limiting situations 
where ip — > 0, we recover the ordinary Bose distribution fk = (e 77 — l) -1 . 
Note that (fk) qP can be rewritten as 

Ufcj9P "g-^i + (i-g)A + ^-gi + (i-^/ fe ( j 

in term of /fc. 

For the purpose of practical applications (in Section 3), it is useful to 
know the development of the distribution (fk) qP in integer series. In this 
respect, we have 

(h) qP = fy^ +i ) ( q — q > + P —A , (6b) 

\q-p p-q J 



3=0 



which reduces to the expansion Ylj=i e m °f the ordinary Bose factor when 
-+ 0. 

3. In the thermodynamical limit, the energy spectrum of the system 
of gp-bosons may be considered as a continuum. Thus, (fk)qp is replaced 
by the ^-dependent factor /(e) which is {fk) qP with = e. Therefore, in 
physical applications, we have to consider integrals of the type 

J s := / e a f(e)de. (11) 
Jo 

By using the development (6b), eq. (11) yields 

J s = T(s + l)(k B T) s+1 a(s + 1) 5P , s > -1, (12) 

where 

oo /Wj'+l) 

+ 1)«p = E fLiw ^' + X ^ " MW- (13) 

j=o u ~ l ~ ^ 



[In eq. (12), T is the Euler integral of the second type.] In the case where 
H = 0, we have a(s + l) qp — > £(s + 1) f° r "~ ** 5 where £(s + 1) is the 
Riemann Zeta series (that converges for s > 0). 

We are now in a position to derive some thermodynamical quantities, 
in D dimensions, for a free gas of N(D) gp-bosons contained in a volume 
V(D). The density p(D) = N(D)/V(D) is given by [14] 



D 
~2 



D 



D 



p(D) = N (D)T(-)(k B T)^a(-) qp , 



(14) 



where 



N (D) := 



1 



D 



2(27r)#r(f + i) g 'n D ' 

In eq. (15), g is the degree of (spin) degeneracy and m the mass of a qp- 
boson. The total energy can then be calculated to be 



m 2 



(15) 



E(D) = N (D)V(D)T(^ + l)(k B T)* +1 v(^ + l)^. (16) 



D 



The specific heat at constant volume easily follows from Cy = (f^)y We 
obtain 



C V (D) = ^N(D)k B 



( 2 +1 > ,(§) 



<7P 



2 7 



9P 



2-(f"l) 



(17) 



Furthermore, the entropy is 



S(D) = N(D)k B 



12 j -(f), 



qp 



k B T 



Note that the state equation for the gas of #p-bosons is 

2 E(D) 

P{D) = dVWy 



(18) 



(19) 



so that the pressure p(D) assumes the same form as in the nondeformed 
case. 

Finally, other thermodynamical quantities may be determined in a sim- 
ple manner by using the thermodynamic potential 

Q(D) = -^N (D)V(D)J» = -^E(D), (20) 

where Jn is given by (12) and (13). 

4. We now examine the condensation of a system of qp-bosons in D 
dimensions. The corresponding density for such a system is given by (14) 
which can be rewritten as 



p(D) =N (D)Jo. 



(21) 



in function of the integral Jd_ v As in the case of ordinary bosons, we 
define the Bose temperature Tb{D) for the gp-bosons by taking \i = 0. By 
introducing \i = in Jd_ 1 , we find that the temperature below which we 
obtain a Bose-Einstein condensation is given by 



where 



T B (D) = ±- 



1 



P(D) 

No(D) r(#) CT0 (f) 



2 /IP J 



_2_ 

D 



D °° i 

<?0(^-)qp = , 77R ([[j + 1]W " Mqp)- 

2 7 =0 + !) 2 



(22) 



(23) 



The study of the convergence of the series (23) determines if Bose-Einstein 
condensation takes place or not. 

In the limiting cases where ip — > 0, eqs. (22) and (23) give back the 
well-known results according to which there is (respectively, there is not) 



condensation in D = 3 (respectively, D = 2) dimensions. As an illustration, 
we recover that the Bose temperature is 



(2tt)^ 2 [ p(3) 



2 
3 



/c^ra |_2.612g 



(24) 



in the case D = 3. 

Returning to the general case, we can prove the following results valid 
for D arbitrary. 

(i) Case (q,p) G R x R : the series (23) converges for < q < 1 
and < p < 1. Thus, there is Bose- Einstein condensation of gp-bosons 
for < q < 1 and < p < 1. In addition, there is no Bose-Einstein 
condensation (i.e., Tb(D) = 0) for g-bosons corresponding to p = 

In contradistinction, Bose-Einstein condensation takes place for g-bosons 
corresponding to p = 1 when < q < 1. 

(ii) Case (q,p) G C x C : Bose-Einstein condensation occurs either for 
the g-bosons corresponding to p = q when < \q\ < 1 or for all g-bosons 
corresponding to p = q = q~ x . 

These results have to be compared to the ones corresponding to the 
limiting situation where (p — > for which Bose-Einstein condensation exists 
only when D > 3. 

5. As a final facet of this Letter, we would like to discuss the con- 
sequence of a gp-deformation of the correlation function of order two 
associated to the radiation field. In the nondeformed case, we know that 
g( 2 ) takes two values, viz., = 1 for a coherent monomode radiation and 
g\ 2 ) = 2 for a chaotic monomode radiation. An interesting question arises : 
Is it possible to interpolate between the latter two values by replacing the 
ordinary bosons of the radiation field by qp-bosons ? 



Our basic hypothesis is to describe the radiation field by an assembly 
of gj)-bosons with the Hamiltonian 



h = ^2h k , h k = Tiu k N k . (25) 
k 

The corresponding Bose statistical distribution for the k mode is then given 
by (10) where rj is replaced by £ = (3TioJk- In the case of gp-bosons, we adopt 
the definition 

(2) = 

{a + a) z 

where a and a + stand for the annihilation and creation operators for the 
k mode. The abbreviation (X) in (26) denotes the mean statistical value 
Z _1 tr (e _/3/l X) for an operator X defined on the Fock space T. Equation 
(26) can be developed as 

9 W =p -i<^>_ (ap) -i<9^. (27 ) 



(a + a) 2 1 ' (a+a) 2 

Of course, (a + a) = (fk) qP as given by (10) with rj = £. In addition, the 
other average values in (27) can be calculated to be 



{[CL a) } ~ (e* - g 2 )(e? - ^)(e£ -p 2 ) 1 * J 



and 



<<a + «> = ? (e ,_; 2)(e ;_ gp) . (29) 



Finally, we obtain 



9 - [q ^ P) eZ-l(et-qi) (e* - «p) (e^ - p 2 ) ' W 

with convergence conditions that parallel the ones for (fk) qP , cf. cases (i)- 
(v) in Section 1 where — must be replaced by £/2. [For instance, when 



{q,p) G R + x R + , we must have < q < es and < p < .] The qp- 
deformed factor depends on the parameters q and p in a symmetrical 
manner (q <-> p symmetry). It also presents a dependence on the energy 
fiujk of the k mode and on the temperature T. In the limiting cases if — > 0, 
we get = 2 that turns out to be the value of the correlation function 
for the choatic monomode radiation of the black body. Let us now examine 
the cases of low temperatures and high energies. In these cases, e^ is the 
dominating term in each of the differences occuring in eq. (30). Therefore, 
we have 

9 {2) ~q + P (31) 

at low temperature or high energy. It is thus possible to reach the value 
g( 2 ) = 1 without employing coherent states. Equation (31) shows that 
we can interpolate in a continuous way from g^ = 1 (coherent phase) to 
g( 2 ) = 2 (chaotic phase). It should be observed that we can even obtain 
either g^ > 2 or g^ < 1 from eq. (30). The situation where g^ < 1 may 
be interesting for describing antibunching effects of the light field (arising 
from the "nonclassical" nature of the radiation field, cf. ref. [21]). 

6. In this work, we concentrated on the influence of a gp-deformation 
of the Weyl-Heisenberg algebra on the statistical properties of a gas of 
bosons in D dimensions. We discussed the phenomenon of Bose-Einstein 
condensation in terms of the deformation parameters q and p. Bose-Einstein 
condensation occurs in most of the domains where the q and p parameters 
may vary. In the case (q,p) G R + x R + with p = q~ x and D arbitrary, 
Bose-Einstein condensation does not take place. For the latter case, we 
may ask the question : Are there any topological interactions (leading to 
some quasi-fermionic behavior through some exclusion rule) that prevent 
the condensation of the corresponding g-bosons ? 



An important result of the present work concerns the correlation func- 
tion g( 2 ) for a gas of photons. The ^-deformation considered here makes it 
possible to obtain a continuum of values for g^ 2 \ In particular, such values 
may be greater than 2 or lesser than 1. 

To close this Letter, it should be observed that alternative choices, non- 
linear in Nk, are possible for the Hamiltonians H (eq. (3)) and h (eq. (25)), 
choices that seem to reflect more deeply the quantum group structure of the 
gp-deformed Weyl-Heisenberg algebra. In this direction, we may mention, 
in the case p = the works of refs. [11,12,15,19,20]. We have justified 
above our choice for eq. (3) (and eq. (4)). The case of eq. (25) is more 
delicate. It would be interesting to investigate the situation where hk = 
(l/2)huk([[Nk\]qp + [[Nk + l]]gp) since the latter expression corresponds to 
a gj)-deformed oscillator. The corresponding partition function can then be 
expressed in series of Bessel functions with some convergence problem. Our 
choice for eq. (25) does not lead to convergence problem and represents 
the first attempt to derive a deformation of the second-order correlation 
function in a nonperturbative way. We hope to return on these matters in 
the future. 
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